Verhulst model with Levy white noise excitation 
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The transient dynamics of the Verhulst model perturbed by arbitrary non-Gaussian 
white noise is investigated. Based on the infinitely divisible distribution of the Levy 
process we study the nonlinear relaxation of the population density for three cases 
of white non-Gaussian noise: (i) shot noise, (ii) noise with a probability density of 
increments expressed in terms of Gamma function, and (iii) Cauchy stable noise. 
We obtain exact results for the probability distribution of the population density 
in all cases, and for Cauchy stable noise the exact expression of the nonlinear relax- 
ation time is derived. Moreover starting from an initial delta function distribution, 
we find a transition induced by the multiplicative Levy noise, from a trimodal prob- 
ability distribution to a bimodal probability distribution in asymptotics. Finally 
we find a nonmonotonic behavior of the nonlinear relaxation time as a function of 
the Cauchy stable noise intensity. 
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I. INTRODUCTION 

The nonlinear stochastic systems with 
noise excitation have attracted extensive at- 
tention and the concept of noise-induced 
transitions has got a wide variety of apphca- 
tions in physics, chemistry, and biology (jH). 
Noise-induced transitions are conventionally 
defined in terms of changes in the number of 
extrema in the probability distribution of a 
system variable and may depend both quan- 
titatively and qualitatively on the character 
of the noise, i.e. on the properties of stochas- 
tic process which describes the noise excita- 
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tion. The Verhulst model, which is a cor- 
nerstone of empirical and theoretical ecol- 
ogy, is one of the classic examples of self- 
organization in many natural and artificial 
systems (0). This model, also known as the 
logistic model, is relevant to a wide range 
of situations including population dynam- 
ics dH; H; 3; H), self-replication of macro- 
molecules (i), spread of viral epidemics (0), 
cancer cell populatiori M) , biological and bio- 
chemical systems (jl; fid ), population of pho- 
tons in a single mode lase r (Illul2 ). autocat- 
alytic chemical reactions (13; 14: 15: [igI: [tT 
freezing of supercooled hquids (jis!), social : 



sci- 



ences (19|; l20t), etc. 



In considering how the population density 
X (t) may change with time t, Verhulst pro- 



posed the following equation 



dx X 
— = rx 1 — — 
dt \ VL 



where there is the Malthus term with the rate 
constant r and a saturation term with the VL 
factor, which is the upper limit for the pop- 
ulation growth due to the availability of the 
resources. 

Really the parameters r and Vl are not 
constant. In fact the parameter r changes 
randomly due to season fluctuations, and the 
parameter VL fluctuates due to the environ- 
mental interaction which causes the random 
availability of resources. As a consequence we 
have the following stochastic Verhulst equa- 
tion 

'" X 



dx 



r (t) X 



1 - 



(2) 



In the context of macromolecular self- 
replication, the model equation ([2]), with con- 
stant VL and a white Gaussian noise in r (t), 
was numerically studied in Ref. (1211 ) and the 
critical slowing down, i.e. a divergence of the 
relaxation time at some noise intensity, was 
found. Later Jackson and co-authors (1221 ) in- 
vestigated the same model, by analog experi- 
ment and digital simulations. They analyzed 
specifically in detail the nonlinear relaxation 
time defined as (1231 ) 



T 



!^[{x{t))-{x{^))]dt 
a;(0) - (a;(oo)) 



(3) 



and did not observe the critical slowing down. 
They explained this discrepancy by the in- 
correct approximate truncation of the asymp- 
totic power series for T used in Ref. (1211). The 



stability conditions were derived in Ref. (|24| ). 
Similar investigations for colored Gaussian 
noise r{t) were performed in Ref. (1251 ) . where 
a monotonic dependence of the relaxation 
time and the correlation time on the noise 
intensity was found. The generalization 
of Eq. ([2]), to study a Bernoulli-Malthus- 
Verhulst model driven by a multiplicative col- 
ored noise, was analyzed recently in Ref. (HT 



The evolution of the mean value in the 
case of Eq. ([2]) with constant r and white 
Gaussian noise excitation BU) = r/fl (t) was 



considered in Refs. (0; M, S; H H; El). In 



Refs. (13; 31) the authors, using perturbation 
technique, obtained the exact expansion in 
power series on noise intensity of all the mo- 
ments and found the long-time decay of t~^/^. 
In Ref. (j3) the authors derived the long-time 
behavior of all the moments of the population 
density by means of an exact asymptotic ex- 
pansion of the time averaged process generat- 
ing function, and found the same asymptotic 
behavior of at the critical point. This 
very slow relaxation of the moments near the 
critical point is the phenomenon of critical 
slowing down. 

In the present paper, using the previously 
obtained results for a generalized Langevin 
equation with a Levy noise source (32; HF 



we investigate the transient dynamics of the 
stochastic Verhulst model with a fluctuating 
growth rate and a constant value for the sat- 
uration population density fl, that is f2 = 
1. The exact results for the mean value of 
the population density and its nonstationary 
probability distribution for different types of 
white non-Gaussian excitation r (t) are ob- 
tained. We find the interesting noise-induced 
transitions for the probability distribution of 
the population density and the relaxation dy- 
namics of its mean value for Cauchy stable 
noise. Finally we obtain a nonmonotonic be- 
havior of the nonlinear relaxation time as a 
function of the Cauchy noise intensity. 



II. STOCHASTIC VERHULST EQUATION 
WITH NON-GAUSSIAN FLUCTUATIONS 
OF GROWTH RATE 

Let us consider Eq. ([2]) with a constant 
saturation value fl = 1, namely 



^ = r (t) X (1 — x) 
dt ^ ^ ^ ' 



(4) 



2 



After changing variable y = ln[a;/(l — x)], we 
obtain 



y{t) = y{0)+ / r{T)dr 
Jo 

and the exact solution of Eq. (jlj) is 
1 - Xo 



xit) 



1 H — exp <,— r (t) dr 

Xo I ^0 



-1 



(5) 

where xq = x(0). Now by substituting in 
Eq. the following expression for the ran- 
dom rate r{t) 



r (t) = r + e (t) 



(6) 



where r > and ^ (t) is an arbitrary white 
non-Gaussian noise with zero mean, we can 
rewrite the solution (151) as 



x{t) 



1 + 



1 - Xo 
( 



-rt-L{t) 



(7) 



Here L [t) denotes the so-called Levy random 
process with L{0) = 0, and £ (t ) = L{t). 



As it was shown in Refs. (I32l : l33l : |3J), Levy 
processes having stationary and statistically 
independent increments on non-overlapping 
time intervals belongs to the class of stochas- 
tic processes with infinitely divisible distri- 
butions. As a consequence, the characteristic 
function of L (t) can be represented in the 
following form (see Eq. (6) in ({3^ ) 




+00 g««^ — 1 — m sin z 



noise. To analyze the time behavior of the 
probability distribution in the transient dy- 
namics it is better not to use the Kolmogorov 
equation for the probability density P{x,t), 
but rather the exact solution ([7]). Using 
the standard theorem of the probability the- 
ory regarding a nonlinear transformation of 
a random variable, we find from Eq. ([7]) 



P{x,t) 



1 



X {1 — x) 





(1 — Xo) X 




-rt,t^ 


Xo (1 - x)_ 







(9) 

where Pl {z, t) is the probability density cor- 
responding to the characteristic function (jS]). 
For a white Gaussian noise £ (t), this distri- 
bution reads 



PL{z,t) 



1 



exp 



-I- 

ADt J 



(10) 



The time evolution of the probability dis- 
tribution P {x,t) for D = 0.5, r = 1, and 
Xo = 0.2 is plotted in Fig. [H 

P{x,t) 




where p {z) is some non-negative kernel func- 
tion. The case p {z) = 2D6 (z) corresponds 
to a white Gaussian noise excitation C, {t), 
while for a symmetric Levy stable noise £ (t) 
with index a we have a power-law kernel 
p{z) = Q with < a< 2. 

In the model under consideration the sta- 
tionary probability distribution has: (i) a sin- 
gularity at the stable point x = 1 for white 
Gaussian noise, and (ii) two singularities at 
both stable points x = and x = 1 for Levy 



(i^fes 1 Time evolution of the probability distri- 
buxion of the population density for white Gaus- 
sian noise excitation with intensity D. The val- 
ues of the parameters are: xq = 0.2, r = 1, 
D = 0.5. 

As it is easily seen, the maximum of the 
unimodal distribution with initial position at 
X = 0.2 shifts with time towards the stable 
point at X = 1. At the same time, as it follows 
from Eqs. ([9]) and (fTOl) . for alH > we have 



lim P (x, t) 



lim P (x, t) 



0. 



:ii) 
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The same picture is observed for an- 
other kernel function p{z) = Kz/ (2sinhz) 
[K > 0), corresponding to a Levy process 
7] {t) with finite moments and the following 
probability density of increments 



PL{z,t) 



/Kt iz 

Tc^r {Kt) VT^7 



y(—-- 

V 2 71 



(12) 

where F (x) is the Gamma function. The cor- 
responding time evolution of the probability 
distribution P (x, t) for K = 0.43, r = 1, and 
Xq = 0.2 is shown in Fig. [2j 



P{x,t) 
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FIG. 2 Time evolution of the probability distri- 
bution of the population density in the case of 
Levy noise with distribution (jl2p . The values of 
the parameters are: xq = 0.2, r = 1, K = 0.43. 

A different situation we have for a Cauchy 
stable noise ^ (t) with constant kernel p (z) = 
Q {a = 1). After evaluation of the integral in 
Eq. (IH]), the probability density of the Levy 
process increments takes the form of the well- 
known Cauchy distribution (js^) 



Dit 



z'' + {D^ty 



(13) 



where = irQ is the noise intensity param- 
eter. In such a case from Eqs. ([9]) and f|T3|) 
for all t > we find 



lim P (x, t) 



lim P 



{x,t) 



oo. 



(14) 



As a result, from an initial delta function 
we immediately obtain a trimodal distribu- 
tion for if: > and then after some transition 



time tc a bimodal one with two singularities 
at the stable points x = and x = 1 (see 
Figs. [3l H and [5]). We should note that the 
transition from trimodal to bimodal distribu- 
tion is a general feature of the model in the 
presence of a Cauchy stable noise, and it is 
\iot limited to some range of parameters. In 
Aact, from Eq. (fT^ and a delta function ini- 
tial distribution inside the interval (0,1), this 
transition always takes place. In the follow- 
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FIG. 3 Time evolution of the probability distri- 
bution of the population density in the case of 
white Cauchy noise excitation. The values of the 
parameters are: xq = 0.2, r = 1, Di = 0.5. 

ing Figs, m and [5] we show the time evolution 
of the probability distribution of the popula- 
tion density for two other values of the noise 
intensity, namely Di = 1 and Di = 5. As the 
noise intensity increases the probability dis- 
tribution shows two singularities near x = 
and X = 1 with different amplitude. 

This transition in the shape of the prob- 
ability distribution of the population den- 
sity is due to both the multiplicative noise 
and the Levy noise source. Using Eqs. ([9]) 
and ffT^ and equating to zero the derivative 
of P{x, t) with respect to x, we obtain the fol- 
lowing condition for the extrema in the range 
< X < 1, and particularly for a minimum 
in the same interval 



z{x, t) 



z{x,t)^ + (D^tf 



X — - 



(15) 
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FIG. 4 Time evolution of the probability distri- 
bution of the population density in the case of 
white Cauchy noise. The values of the parame- 
ters are xq = 0.2, r = 1, Di = 1. 
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FIG. 5 Time evolution of the probability distri- 
bution of the population density in the case of 
white Cauchy noise. The values of the parame- 
ters are xq = 0.2, r = 1, Di = 5. 



with 



z{x, t) = In 



Xo)x 



xo (1 - x) 



rt . 



(16) 



This condition can be solved graphically by 
finding the intersection between the func- 
tions yi = z{x,t)/{z{x,tY + {DitY) and 
y2 = X — 1/2. This is done in the follow- 
ing Figs. El El where the function yi is 
plotted for three different values of time and 
noise intensity. In each figure the black blue 
curve (color on line) corresponds to the criti- 
cal value of time tc for which we have a noise 
induced transition of the probability distribu- 
tion of the population density from trimodal 



to bimodal, that is from two minima and one 
maximum to one minimum inside the interval 
< a; < 1. The appearance of one minimum 
in the probability distribution is the signa- 
ture of this transition. The three values 




FIG. 6 Plots of both sides of Eq. (white 
Cauchy noise): function yi (solid curves), func- 
tion 2/2 (dashed curve) , for three values of time, 
namely: t = 1, 1.5, 2.493. The critical time is 
tc = 2.493 (black blue curve). The values of the 
other parameters are: xq = 0.2, r = 1, Di = 0.5. 




FIG. 7 Plots of both sides of Eq. ([T5]) (white 
Cauchy noise): function yi (solid curves), func- 
tion 2/2 (dashed curve) , for three values of time, 
namely: t = 0.5, 1, 1.635. The critical time is 
tc = 1.635 (black blue curve). The values of the 
other parameters are: xq = 0.2, r = 1, Di = 1. 

of the critical time tc corresponding to the 
three values of the Levy noise intensity in- 
vestigated are: Di = 0.5, tc = 2.493; Di = 
1, tc = 1.635; Di = 5,tc = 0.247. One rough 
evaluation of the critical time tc is obtained 
by putting equal to 1 the scale parameter of 
the Cauchy distribution of Eq. (fT3l) . that is 



5 



tc ~ I/-D1. The critical time tc is the time at 
which the maximum and one minimum of the 
probabihty distribution (see Figs. [3], HI and [5]) 
coalesce in one inflection point and in this 
point X the function 7/2 = x — 1/2 becomes 
tangent at the function yi (see Figs. O El 
and [8]) . It is interesting to note that the crit- 
ical time tc decreases with the noise intensity 
Di. This is because by increasing the noise 
intensity, more quickly the population den- 
sity reaches the two points near the bound- 
aries X = and x = 1. 




FIG. 8 Plots of both sides of Eq. ([I5]) (white 
Cauchy noise): function yi (solid curves), func- 
tion 2/2 (dashed curve), for three values of time, 
namely: t = 0.5, 1, 1.635. The critical time is 
tc = 0.247 (black blue curve). The values of the 
other parameters are: xq = 0.2, r = 1, Di = 5. 



III. NONLINEAR RELAXATION TIME OF 
THE MEAN POPULATION DENSITY 

It must be emphasized that to find the 
time evolution of the mean population den- 
sity one can use two different approaches. 
The first one was proposed in Ref. (1221 ) . Ac- 
cording to the exact solution ([7]) of the Ver- 
hulst equation (jlj), we can rewrite this ex- 
pression in the following form 



x{t) = f (e-'-*-^^*)) 



(17) 



where 



1 - Sn ^ -1 



Xo 



Then, by expanding the smooth function f|T8|) 
in a standard Taylor power series in q around 
the point g = we have 



n=0 



(19) 



After substitution of Eq. ^ in Eq. ^ and 
averaging we obtain 



(X 



(t)) = E 



^ (0)e- 



■nrt 



n=0 



-nL{t) 



(20) 



or, in accordance with Eq. ([8]) 



n=0 "'■ 



■nrt 



X exp It J 



-co Q _ X _|_ gin z 



(21) 



p (z) dz > . 



For white Gaussian noise ^ (t) with kernel 
p (z) = 2D6{z) we obtain from Eq. fl2ip the 
following asymptotic series 

W) = E f^e^*" (22) 



n=0 



By considering a finite number of terms in 
this expansion leads to a wrong conclusion 
about the critical slowing down phenomenon 
in such a system, as found in Ref. (l2ll ). The 
exact result is obtained, of course, by sum- 
ming all the terms in Eq. fl22|) . Moreover, 
for most of the kernels p{z) the integral in 
Eq. fl2Tl) diverges. Thus, this approach is in- 
appropriate for our purposes, and it is better 
to use the direct average in Eq. (iTl). There- 
fore, using this second approach we have 



(x(t)) 



1 + l^e-*-^ 

Xo 



Pl {z, t) dz. 
(23) 



Let us consider now different models of 
white non-Gaussian noise ^ (t). We start with 
the white shot noise 



(24) 
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having the symmetric dichotomous distri- 
bution of the pulse amphtude -P(a) = 
[S {a — ao) + 6 {a + oq)] /2, mean frequency z/ 



of pulse train, and kernel p (z) 
From Eq. (jSj) we have 



-ut{l—cos aou) 



uz'P iz) 



(25) 



By making the reverse Fourier transform in 
Eq. fl25|) we find the probability distribution 
of the corresponding Levy process 



naQ) 



(26) 

where /„ (x) is the n-order modified Bessel 
function of the first kind. The relaxation of 
the mean population density x{t) is shown in 
Fig.H According to the Eq. and ([26]) the 
stationary value of the population density in 
such a case is {x)^^ = 1, but the relaxation 
time ([3]) increases with increasing the mean 
frequency of pulses. 



<x(0> 




t 

7 

FIG. 9 Nonlinear relaxation of the mean popu- 
lation density in the case of white shot noise ex- 
citation, for three values of the mean frequency 
v, namely i' = 0.1, 1,3. The values of the other 
parameters are: xq = 0.2, r = l,ao = 1- 

For white non-Gaussian noise with the 
kernel p{z) = Kz/ (2sinh2;) we observe a 
similar transient dynamics, which is shown in 
Fig. [TOl We have the same stationary value 
and the relaxation time T increases 



<x(t)> 




t 

1 

FIG. 10 Nonlinear relaxation of the mean pop- 
ulation density in the case of Levy noise with 
distribution ()13p , for three values of the param- 
eter namely K = 0.2, 1, 5. The values of the 
other parameters are: xq = 0.2, r = 1. 

with increasing the parameter K, which is 
proportional to the noise intensity. 

Finally, in the case of white Cauchy noise 
^ {t) we obtain interesting exact analytical re- 
sults. First of all, substituting Eq. ( fT3l) in 
Eq. fl2^ and changing the variable z = Dity 
under the integral, we obtain 



(x(t)) 



1 



+ 00 ' 



1 



Xq 



dy 



{x) 



St 



(27) 
we find 



For the stationary mean value 
from Eq. ([27D 

1 (r + Diy) dy 



{x) 



St 



1 

lim (x (t)) = - / 



(28) 

where 1 (x) is the step function. After eval- 
uation of the integral in Eq. (!28|) we obtain 
finally 



St 



11 r 
— I — arctan — . (29) 

As it is seen from Fig. [TT] and Eq. ( !29l) . 
for small noise intensity Di, with respect 
to the value of the rate parameter r = 1, 
the stationary mean value of the population 
density is approximately 1, as for the other 
white non-Gaussian noise excitations consid- 
ered. But for large values of Di, this asymp- 
totic value, which is independent from the 
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initial value of population density xq, tends 
to 0.5. 



T 



1 u ^i^^r" 




t 
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FIG. 11 Nonlinear relaxation of the mean popu- 
lation density in the case of white Cauchy noise, 
for four values of the noise intensity -Di, namely 
Di = 0.1,0.5,1,5. The values of the other pa- 
rameters are: xq = 0.2, r = 1. 

It is interesting also to analyze, for this 
case of white Cauchy noise, the dependence 
of the relaxation time T from the noise inten- 
sity Di. Substituting Eq. ([27]) in Eq. ([3]) and 
changing the order of integration, for initial 
condition xq = 0.5, we are able to calculate 
analytically the double integral in t and in y 
obtaining the final result 



r {1 + Dj/r'^) arccot{Di/r) 

We find a non-monotonic behavior of the re- 
laxation time T versus the noise intensity Di 
with a maximum at the noise intensity Di = 
0.43, as shown in Fig. [121 This nonmonotonic 
behavior is also visible for another initial po- 
sition xo = 0.2, as shown in Fig. [TTl Here the 
relaxation time to reach the stationary value 
of population density Xst increases from very 
low noise intensity {Di = 0.1) to moderate 
low intensity {Di = 0.5), while decreases for 
higher noise intensities {Di = 1.5). This is 
also due to the dependence of Xst from the 
noise intensity Di (see Eq. ( [30]) ). We note 
that this non-monotonic behavior of the re- 
laxation time T is related to the peculiarities 



1.5 
1.25 
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FIG. 12 Nonmonotonic behavior of the nonlin- 
ear relaxation time T as a function of the white 
Cauchy noise intensity Di. The values of the 
other parameters are: xq = 0.5, r = 1. 

of the transient dynamics of the mean popu- 
lation density and it will be object of further 
investigations. 

IV. CONCLUSIONS 

The transient dynamics of the Verhulst 
model, perturbed by arbitrary non-Gaussian 
white noise, is investigated. This well-known 
equation is an appropriate ecological and bi- 
ological model to describe closed-population 
dynamics, self-replication of macromolecules 
under constraint, cancer growth, spread of vi- 
ral epidemics, etc... By using the properties 
of the infinitely divisible distribution of the 
generalized Wiener process, we analyzed the 
effect of different non-Gaussian white sources 
on the nonlinear relaxation of the mean pop- 
ulation density and on the time evolution of 
the probability distribution of the popula- 
tion density. We obtain exact results for the 
nonstationary probability distribution in all 
cases investigated and for the Cauchy stable 
noise we derive the exact analytical expres- 
sion of the nonlinear relaxation time. Due to 
the presence of a Levy multiplicative noise, 
the probability distribution of the popula- 
tion density exhibits a transition from a tri- 
modal to a bimodal distribution in asymp- 
totics. This transition, characterized by the 
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appearance of a minimum, happens at a crit- 
ical time tc, which can be roughly evaluated 
as tc ~ (where Di is the noise inten- 

sity) and exactly evaluated from the condi- 
tion f[T^ . Finally a nonmonotonic behavior 
of the nonlinear relaxation time of the pop- 
ulation density function of the Cauchy 
noise intensity was found. 
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